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Abstract 

A lattice Delaunay polytope D is called perfect if it has the prop- 
erty that there is a unique circumscribing ellipsoid with interior free 
of lattice points, and with the surface containing only those lattice 
points that are the vertices of D. An inhomogeneous quadratic form 
is called perfect if it is determined by such a circumscribing "empty 
ellipsoid" uniquely up to a scale factor. Perfect inhomogeneous forms 
are associated with perfect Delaunay polytopes in much the way that 
perfect homogeneous forms are associated with perfect point lattices. 
We have been able to construct some infinite sequences of perfect De- 
launay polytopes, one perfect polytope in each successive dimension 
starting at some initial dimension; we have been able to construct an 
infinite number of such infinite sequences. Perfect Delaunay poly- 
topes are intimately related to the theory of Delaunay polytopes, and 
to Voronoi's theory of lattice types. 

1 Introduction 

Consider the lattice Z*^, and a lattice polytope D. If D can be circumscribed 
by an ellipsoid £ = c}{ x G M"^ | /(x) < -R^}, where / is a quadratic function 
in with no Z'^-elements interior to { x G vert /(x) < -R^}, 

and with all Z'^-elements on £ being the vertices of -D, we will say that D 
is a Delaunay polytope with respect to the metric form defined by £; more 
informally, we will say that D is Delaunay if it can be circumscribed by 
an empty ellipsoid £. Typically there is a family of empty ellipsoids that 
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can be circumscribed about a given Delaunay polytope D but, if there is 
only one, so that £ is uniquely determined by D, we will say that D is a 
perfect Delaunay polytope. Perfect Delaunay polytopes are distinguished in 
that they cannot fit properly inside other lattice Delaunay polytopes. Perfect 
Delaunay polytopes are fascinating geometrical objects - examples are the 
six and seven dimensional Gossett polytopes with 27 and 56 vertices known 
as 221 and 821 in Coxeter's notation, which appear in the Delaunay tilings 
for the quadratic forms corresponding to the root lattices and Ey. 

We have studied perfect Delaunay polytopes by constructing infinite se- 
quences of them, one perfect Delaunay polytope in each successive dimension, 
starting at some initial dimension. We have been able to construct an infi- 
nite number of infinite sequences of perfect Delaunay polytopes. One of our 
constructions is the sequence of G-topes, G^, d = 6,7,... , with the initial 
term being the six- dimensional Gossett polytope 221 with 27 vertices; each 
G-tope is asymmetric with respect to central inversion, and G'' has ("^2^) — 1 
vertices. Another construction is the sequence of C-topes, C^, d — 7,S, ... , 
with initial term the seven-dimensional Gossett polytope 831 with 56 vertices; 
each C-tope is symmetric with respect to central inversion, and has 2 ('^^^) 
vertices. Just as the six-dimensional Gossett polytope can be represented as 
a section of the seven dimensional one, each term G*^ of the asymmetric se- 
quence can be represented as a section of the term G*^"*"^ of the symmetric 
sequence. 

Each perfect Delaunay polytope in the sequence of G-topes uniquely de- 
termines a lattice that is an analogue of the root lattice Eq, and each term 
in the sequence of G-topes uniquely determines a lattice that is an analogue 
of the root lattice £'7. The Voronoi and Delaunay tilings for the lattices in 
these sequences of lattices show many features of the lead terms, namely, the 
Voronoi and Delaunay tilings for the quadratic forms corresponding to the 
root lattices Eq and E^. 

There is a number of reasons why perfect Delaunay polytopes are fas- 
cinating objects for study. First, as mentioned before, they are "perfect" 
inhomogeneous analogs of perfect lattices. This alone seems to be a natural 
motivation for a geometer of numbers. If the empty ellipsoid £ circumscrib- 
ing a perfect Delaunay polytope D with the vertex set vert D is defined by 
an equation /(x) = 0, where / is an inhomogeneous quadratic form, then 
all the coefficients of / are uniquely determined (up to a common scaling 
factor) from the system {/(x) = 1 | x e vertD}. Glearly, vertD is then the 
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set of all points on which / attains its minimum over Z'^. An analogous prop- 
erty for homogeneous quadratic forms is called perfection after Korkin and 
Zolotareff (1873) (see Martinet (2003) and Conway and Sloane (1988) for 
modern treatment of perfect homogeneous forms). Naturally, in our context 
/ is called an inhomogeneous perfect form and £ is referred to as a perfect 
ellipsoid. The vertices of a perfect Delaunay polytope are analogs of the 
minimal vectors for a perfect form: the minimal possible number of vertices 
of a perfect Delaunay polytope is "^^"^^^^ -)_ while the minimal number of 
shortest vectors of a perfect forms is !ll!Ltl)_ 

Perfect Delaunay polytopes were first considered by Erdahl in 1975 in con- 
nection with lattice polytopes arising from the quantum mechanics of many 
electrons. He showed (1975) that there are perfect Delaunay polytopes in 
one-dimensional lattices, and showed that the Cosset polytope C^ = 22i 
with 27 vertices was a perfect Delaunay polytope in the root lattice E^. He 
also showed that there were no perfect Delaunay polytopes in dimensions 2, 
3, and 4. These results were extended by Erdahl (1992) by showing that 
the 7-dimensional Cosset polytope C^ = 331 with 56 vertices is perfect, and 
that there are no perfect Delaunay polytopes in lattices with dimension lying 
between one and six. Erdahl also proved that [0, 1], C^, and C'' are the only 
perfect Delaunay polytopes existing in root lattices. Deza, Crishukhin, and 
Laurent (1992,1995) found more examples of perfect Delaunay polytopes in 
dimensions 15, 16, 22, and 23. The first construction of infinite sequences 
of perfect Delaunay polytopes was first announced at the Conference ded- 
icated to the Seventieth Birthday of Sergei Ryshkov (Erdahl, 2001), and 
later reported by Rybnikov (2001) and Erdahl and Rybnikov (2002). Perfect 
Delaunay polytopes have been classified up to dimension 7 - Dutour (2004) 
proved that C® = 221 is the only Cosset polytope for d = Q. It is strongly sus- 
pected that the existing lists of seven and eight dimensional perfect Delaunay 
polytopes are complete (see http://www.liga.ens.fr/~dutour/). 

Perfect Delaunay polytopes are important in the study of strongly regular 
graphs or, more generally, theory of association schemes. In fact, the Schlafli 
graph and the Cosset graph can be constructed as 1-skeletons of Cosset 
polytopes C^ and C'', and that is precisely how the Cosset graph has been 
discovered. A maximal family of equiangular lines is a metrical concept, 
closely related to the combinatorial notion of strongly regular graph (see 
Deza et al (1992), Deza and Laurent (1997) for details). All Delaunay 
polytopes found by Deza et al have been derived from such families of lines. 
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Perfect Delaunay polytopes play an important role in the L-type reduc- 
tion theory of Voronoi and Delaunay for positive (homogeneous) quadratic 
forms. An L-type domain is the collection of all possible quadratic metric 
forms that give the same Delaunay tiling D for T/'. L-type domains are 
relatively open polyhedral cones, with boundary cells that are also L-type 
domains - these conical cells fit together to tile the cone of metrical quadratic 
forms, which is described in the next section in detail. Simplicial Delaunay 
tilings label the full dimensional conic "tiles", and all other possible Delau- 
nay tilings label the lower dimensional cones. A significant role for perfect 
Delaunay polytopes is that they provide labels for a subclass of edge forms 
for L-types. Edges-forms are forms lying on extreme rays of full-dimensional 
L-type domains, which are labels by simplicial Delaunay tilings. Prior to the 
discovery of the sequence {G'^} of G-topes by Erdahl and Rybnikov in 2001 
only finitely many edge forms were known. All of the infinite sequences 
of perfect Delaunay polytopes that we report on correspond to infinite se- 
quences of edge forms for L-type domains. The significance of extreme L- 
types is much due to their relation to the structure of Delaunay and Voronoi 
tilings for lattices. The Delaunay tilings that correspond to edge forms play 
an important role: The Delaunay tiling for an L-type domain is the inter- 
section of the Delaunay tilings for edge forms for the L-type. (Erdahl, 2000) 
There is a corresponding dual statement on the structure of Voronoi poly- 
topes: The Voronoi polytope for a form ip contained in an L-type domain 
is a weighted Minkowski sum of linear transforms of Voronoi polytopes for 
each of the edge forms. The latter dual result was first established by H.-F. 
Loesch in his 1990 doctoral dissertation, although it was first pubhshed by 
Ryshkov (1998, 1999), who independently rediscovered Loesh's theorem; this 
dual result was given a shorter and simpler proof by Erdahl (2000). 

Edge forms that are interior to the cone of metric forms are rare in low 
dimensions. They first occur in dimension 4: has an extreme L-type. A 
good proportion, but not all, of the edge forms appearing in lower dimensions 
relate either directly or indirectly to perfect Delaunay polytopes. As shown 
by Dutour and Vallentin (2003) this situation does not persist in higher 
dimensions: there is an explosion of edge forms in six dimensions, and that 
only a tiny fraction of these are inherited from perfect Delaunay polytopes. 

All of our sequences of perfect Delaunay polytopes determine correspond- 
ing sequences of lattices with similar combinatorial properties. Our construc- 
tions are first steps of a program to explore the geometry of higher dimen- 
sional lattices through infinite sequences of lattices. The infinite sequences 
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we have constructed are particularly interesting because of the role they play 
in the structure theory of Delaunay polytopes and Delaunay tilings, which 
are described in the following section. 

2 Homogeneous and Inhomogeneous 
Forms on Lattices 

The Voronoi and Delaunay tilings for point lattices are constructed using the 
Euclidean metric, but are most effectively studied by mapping the lattice 
onto Z*^, and replacing the Euclidean metric by an equivalent metrical form. 
For a d- dimensional point lattice A with a basis bi, b2, this is done as 
follows. If V is a lattice vector with coordinates Zi, Z2, Zd relative to this 
basis, then v can be written as v = Bz, where B = [bi, b2, bd] is the basis 
matrix. The squared Euchdean length, is given by |vp = z^B^Bz =(^b(z), 
where z is the column vector given by [zi, Z2, Zd]'^ . This squared length 
can equally well be interpreted as the squared length of the integer vector 
z eZ'^ relative to the metrical form ip-^. Therefore, the Voronoi and Delaunay 
tilings for A, constructed using the Euclidean metric, can be studied using 
the corresponding Voronoi and Delaunay tihngs for Z'' constructed using the 
metrical form (p-^. Moreover, variation of the Voronoi and Delaunay tilings 
for A in response to variation of the lattice basis bi, b2, b^ can be studied 
by varying the metrical form ip for the fixed lattice Z'^. In the discussion 
below we will keep the lattice fixed at Z'^, and vary the metrical form (p. 
With slight abuse of terminology we call any semidefinite form metric. 

Definition: The squared length of a vector v with respect to a metric form 
(p is called the norm ofv relative to (p. 

The inhomogeneous domain of a Delaunay polytope: The following 
discussion of the geometry of J" is a summary of some results contained in 
(Erdahl, 1992). Let T be the cone of quadratic functions on M'^ defined by: 

? = {/ e R[xi, . . . I deg/ = 2, Vz e /(z) > 0}. 

The condition Vz e /(z) > requires the quadratic part of / to be 
positive semi-definite, and requires any subset of R'^ where / assumes negative 
values to be free of Z'^-elements. The real surface determined by the equation 
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/(x) = might be empty set, it might be a subspace, or it might have the 
form 

S.f = EqX K, 

where Sq is an elhpsoid and K a complementary subspace. The last case is 
the interesting one - depending on the dimension oi K, £/ is either an empty 
(of lattice points) ellipsoid or an empty cylinder with ellipsoid base. 

We denote by V{f) the set £/ n Z'^. In the case where the surface 
includes integer points and is an empty ellipsoid, V{f) — Ef nZ'^ is the vertex 
set for the corresponding Delaunay polytope Df = conv\/(/). Convcrsly, 
if D is a Delaunay polytope in Z"', there is a circumscribing empty ellipsoid 
8,D determined by a function G CP. More precisely, since D is assumed 
Delaunay there is a metrical form c/?^, a center c and radius R, so that 
/d(x) = (pjj{x — c) — — is the equation of a circumscribing empty 
ellipsoid Ed- Since //) is non-negative on Z"', it is an element of T. 

In the case where £/ is an empty cylinder, there will be an infinite number 
of integer points lying on this surface. When this happens V{f) = £/ fl Z'^ 
is the vertex set for a non-bounded Delaunay polyhedron Df — couvV{f), 
which is a cylinder with a base, which is also called a Delaunay polytope. 
Recall than a convex polyhedron is called a polytope when it is bounded. 

Definition : Let D be a Delaunay polyhedron. Then, the (inhomogeneous) 
domain for D is: 

= T I Df^D} 

Such domains are relatively open convex cones that partition int T (in fact, 
they partition a larger subset of J", which includes int "?) . 

The elements / e satisfy the homogeneous linear equations /(z) = 0, z e 
vert D — D fl Z*^, so the dimensions of domains vary depending upon the 
rank of this system. When D is a single edge, the rank is one and I'd is a 
relatively open facet of the partition with dimension dim CP — 1 = ('^^^) — 1. 
When the rank is equal to ('^^^) — 1, which is the maximum possible in order 
that To not be empty, T/j is an extreme ray of CP. 

Definition: Let V{p) be the set of integer points lying on the boundary 
of a d-dimensional Delaunay polyhedron. A function p e CP is perfect if and 
only if the equations p(z) = 0, z & V(p) uniquely determine p up to 
scaling. In this case we will call the subset V(j)) C Z*^ perfect, and will call 
Dp — conv V{p) perfect. 



6 



The elements of a 1-dimensional inhomogeneous domain are perfect, and the 
Delaunay polytopes that determine such domains are perfect. The perfect 
subsets V{p) must be maximal among the subsets { V{f) C | / G CP}. 

Perfect inhomogeneous forms are analogues of perfect homogeneous forms — 
both achieve their minimum value on Z'^ sufficiently often that the represen- 
tations of the minimum determine the form. In the homogeneous case the 
minimum is taken on the non-zero elements of Z*^; if the minimum value 
is achieved on the set V{lp), then (/? is uniquely determined by the equations 
(p{z) = m^p, z ^V{(p). Similarly, in the inhomogeneous case the minimum 
is taken on Z*^, and the minimum value is zero; if this minimum value is 
achieved on the set V{f) C Z'^, then / is uniquely determined (up to a scale 
factor) by the equations /(z) = 0, z eV{f). 

Theorem 1 If p E is perfect, then there is an ellipsoid 8,q, where < 
dim(£o) < d, and complementary subspace K so that Ep — Eq x K, and so 
that £o and K satisfy the following arithmetic conditions: 

£o n Z*^ are the vertices of a perfect Delaunay polytope in (aff £o) ^ Z*^; 
Kr\Z'^ is a sublattice such that Z'^ = ((aff £o) n Z'^) © (Xn Z*^). 

In this case we have V{p) = £p fl Z"' = (£o fl 71^) 0(A' fl Z^). Conversely, 
any ellipsoid £o and an affine complementary subspace K , satisfying these 
arithmetic conditions, determine a surface 8.0 x K for a perfect element p of 
y and, therefore, determine a perfect element up to a scale factor. 

By this theorem, if p is perfect, then is either a perfect Delaunay poly- 
tope or, in the degenerate case when dimX > 0, a cylinder with a perfect 
Delaunay polytope as base. 

A Z'^-vector [-Ui, . . . is called primitive if GCD(-Ui, . . . = 1. As 
an example, consider primitive vectors a, b G Z'^ such that a • b =1. Then 
£0 = {0, b} is an empty ellipsoid in the 1-dimensional lattice (aff £0) n Z*^, 
and K — a.^ — {xeR'^|x-a = 0} is a complementary subspace. The 
subset {0, b} = £0 n Z'^ is also the vertex set for a Delaunay polytope in 
aff (£onZ''), and is defined so that Z'^ = ((aff £0) nZ'')) 0(/s:nZ"'). These 
are the arithmetic conditions stated in the theorem, so that {0, b} x K is 
the surface for a perfect element p of IP, which, up to a scale factor, is given 
by p(x) = (a • x)(a • x— 1). The preimages of negative real values of p lie 
between two hyperplanes with equations a • x =0 and a • x =1, a region which 
is a degenerate Delaunay polyhedron. 
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With the exception of the one-dimensional perfect domains, all inhomo- 
geneous domains CP/j have proper faces that are inhomogeneous domains of 
lesser dimensions. If D is inscribed into another Z'^-polytope D' so that 
verti:'' D verti:', then J^^' C d'^o- 

In the following definition the Delaunay polyhedron D may be of any 
dimension. 

Definition: Let D be a lattice Delaunay polyhedron. Then, the inhomo- 
geneous domain for D is defined as: 

Tz? = {/eT I vertD = 



If D is a d-dimensional Delaunay simplex then dim?/) = ('^g^), but if D is 
a perfect Delaunay polyhedron, then dimlPc = 1. 

Theorem 2 Let D he a d- dimensional Delaunay polyhedron. Then 

= { iopp iOp e M>o }, 

where the summation is over all perfect elements p such that vert D C V{p). 

In general, not all relatively open faces on the boundary of an inhomo- 
geneous domain CP/j are inhomogeneous domains. However, in the special 
case where D is d- dimensional, all extreme rays are perfect inhomogeneous 
domains and all relatively open faces are inhomogeneous domains (see Er- 
dahl 1992). In this case an arbitrary element / G To has the following 
representation: 

/ = ^ UpP, 

{P I V(p)5vertD} 

where cjp > 0. The summation is over the perfect elements p with the 
property that vert D C V{p). 

These last two paragraphs, and our representation theorem for Delaunay 
polytopes, show the important role played by perfect Delaunay polytopes in 
the theory. 

The homogeneous domain of a Delaunay tiUng: Voronoi's classifica- 
tion theory for lattices, his theory of lattice types (L-types), was formulated 
using metrical forms and the fixed lattice Z"'. In this theory two lattices are 
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considered to be the same type if their Delaunay tihngs are affinely equiva- 
lent. Consider a positive definite quadratic form kp. Then a lattice polytope 
D is Delaunay relative to (/? if it can be circumscribed by so called empty 
ellipsoid £ with equation of the form 

</7(x — c) = i?^, 

where c e M*^ and R G M>o; in addition, conv £ must have no interior 
Z'^-elements, and verti? must be given by £nZ'^. The collection of all 
such Delaunay polytopes fit together facet-to-facet to tile W^^ a tiling that is 
uniquely determined by ip. This is the Delaunay tiling for relative 
to the metrical form ip. If a second metrical form d has Delaunay tiling 
Dij, and if T)^ is identicle to, or affinely equivalent to Di^, then ip and are 
metrical forms of the same L-type. 

The description we give below requires that certain degenerate metrical 
forms be admitted into the discussion, namely, those forms ip for which K = 
Keiip is a rational subspacc of W^. The Delaunay polyhcdra for such a 
form are themselves degenerate-they are cylinders with axis K and Delaunay 
polytopes as bases. These cylinders fit together to form the degenerate 
Delaunay tiling T)^. For example, if a eZ'' is primitive, </^(x) = (a • x)^ is 
such a form — the kernel K is the solution set for ip{'x) =0, and given by a-*-, 
which is rational. The Delaunay tiles arc infinite slabs, each bounded by a 
pair of hyperplanes a ■ x =k, a ■ x =k + 1, k E Z. These fit together to tile 

Let $ be the cone of metrical forms in d variables, namely, the cone 
of positive definite quadratic forms and semidefinite quadratic forms with 
rational kernels. For each metrical form e $ there is a Delaunay tiling 

for Z'^. 

Definition: If D is a Delaunay tiling for Z"^, then the following cone of 
positive definite quadratic forms 

$D = { e I = D } 

is called an L-type domain. 

For this definition the Delaunay tilings can be the usual ones, where the tiles 
are Delaunay polytopes - or they could be degenerate Delaunay tilings where 
the tiles are cylinders with a common axis K. 
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The relatively open faces of an L-type domain, are L-type domains. If 
D is a trangulation, $xi has full dimension ('^^^) ; this is the generic case. If 
D is not a triangulation, dim$D is less than ('^^^), and $xi is a boundary 
cell of a full-dimensional L-type domain. In the case where dim$D = 1, the 
elements </? e $d are called edge forms as they correspond to extreme rays 
of full-dimensional L-type domains. 

Let D be a polytopc in the Dclaunay tiling D. Then, if 7r$ is the 
projection onto the quadratic part, $d C 7iq,{7u). Since this containment 
holds for all Delaunay tiles D E D, there is the following description of $d 
in terms of inhomogeneous domains: 

where the intersection is over all (i-dimensional Delaunay polytopes in T). 
Since the containment $d C 7:^(7 d) also holds for all Delaunay tilings D 
that contain D, there is the following description of 7r$(!P£)) in terms of 
homogeneous domains: 

where the disjoint union is over all Delaunay tilings for Z'^ that contain 
D. This holds not only for full-dimensional cells oi Td, biit for cells of all 
dimensions. The last equality shows that 7:^(^0) is tiled by L-type domains. 
It also establishes the following 

Theorem 3 Ifp E 7 is perfect and, therefore, an edge form for an inhomo- 
geneous domain, then 7r$(p) is an edge form for an L-type domain. 

This result can also be established in a more direct way by appealing to the 
definition of L-type domain: if D is a perfect Delaunay polytope, or perfect 
degenerate Delaunay polyhedron, then 7r$(yz)) is a one-dimensional L-type 
domain. By definition, the elements of 7r$(3'D) are then edge forms. 

The converse of theorem does not hold - there are edge forms for L-type 
domains that are not inherited from perfect elements in CP. Evidence has 
accumulated recently that the growth of numbers of types of edge forms with 
dimension is very rapid starting in six dimensions (see Dutour and Vallentin, 
2003), but the the growth of perfect inhomogeneous forms is much less rapid 
- there is some hope that a complete classification can be made through 
dimension nine, or even through ten. 
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3 Infinite Sequences of 

Perfect Delaunay polytopes 

Consider the following sets of vectors in R*^: 

Diu = { [l^ ^'-'] - } u { [r^\ Q"—'] - } 

where s,k E N and j — [1^] (1^ means the entry 1 is repeated d times, and 
similarly for 1* and O*^"*). All permutations of entries are taken so that 
= (j) + (^^J = ('^+1). The following is the main theorem of this 

paper. 

Theorem 4 For d > k{2s + 1) + 1, where s > 1, k >2, the polytope 

nS = ^conv{ L>J, U -L>f,, }. 

is a symmetric perfect Delaunay polytope for the affine lattice Af ^, = aff 
the origin is the center of symmetry for P^j^ and does not belong to Af^^,. The 
circumscribing empty ellipsoid can be defined by the equation ^.(x) = R^, 
where 

(p^ki^) = 4k{d - k{2s + l))|x|2 + ((i^ - {4k + 2s+ l)d + 4k{2s + A;))(j • x)l 

Each pair of positive integers (s, A;), for s > 1, A; > 2, determines an infinite 
sequence of symmetric perfect Delaunay polytopes, one in each dimension, 
with the initial dimension given by k{2s + 1) + 1. For s = 1, k = 2 the infinite 
sequence is the one described in the opening commentary, i.e., C°',(i > 7, 
where the initial term is the Gosset polytope 821 = with 56 vertices. 

The ('^^^) diagonal vectors D'^j^ for P^ ,^ have the origin as a common mid- 
point, forming a segment arrangement that generalizes the cross formed by 
the diagonals of a cross-polytope. Moreover, these {^~^^) diagonals are prim- 
itive and belong to the same parity class for Af namely, they arc equivalent 
modulo 2A^^^. More generally, primitive lattice vectors u, v in some lattice 
A, with mid-points equivalent modulo A, are necessarily equivalent modulo 
2A and thus belong to the same parity class. And conversely, the mid-points 
of lattice vectors u, v G A belonging to the same parity class are equivalent 
modulo A. By analogy with the case of cross-polytopes, we call any such 
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arrangement of segments or vectors a cross. The convex hulls of such crosses 

often appear as cells in Dclaunay tilings - cross polytopes are examples, as 
are the more spectacular symmetric perfect Delaunay polytopes. There is 
a criterion due to Voronoi (1908) and Baranovskii (1991) that determines 
whether these crosses are Delaunay: Let A be a lattice, let ip be a metrical 
form, and let C be the convex hull of a cross of primitive vectors belonging 
to the sam,e parity class. Then C is Delaunay relative to ip if and only if 
the set of diagonal vectors forming the cross is the complete set of vectors of 
minimal length, relative to ip, for their parity class. This is the criterion 
we have used to establish the Delaunay property for the symmetric perfect 
Delaunay polytopes P^^. 

The following result shows that asymmetric perfect Delaunay polytopes 
can appear as sections of symmetric ones. 

Theorem 5 Ford>6 letu^ [-1^, I'^-i] e Z'^+^ Then 

= conv{ V e vert P^^^ \ v • u = ^ } 

is an asymmetric perfect Delaunay polytope for the affine suhlattice Af ^ = 
affvertG"^. The circumscribing empty ellipsoid can be defined as {x e 
affG^ I (^^J^(x) = R^}, where 

(^f+i(x) = 8((i - 5)|x|2 + (^2 _ 9(i + 22)(j • x)'. 

This is the infinite sequence of G-topes described in the introduction, with 
Gossett polytope 22i = G^ as the initial term for d — 

The terms in this sequence have similar combinatorial properties. For 
example, the lattice vectors running between vertices all lie on the boundary, 
in all cases. These lattice vectors are either edges of simplicial facets, or 
diagonals of cross polytope facets-there are two types of facets, simplexes and 
cross polytopes. The 6-dimensional Gossett polytope has 27 five-dimensional 
cross-polytopal facets, but for the rest the number is twice the dimension, 
2d. can be found as a section of G^, but G® and G^ do not have sections 
arithmetically equivalent to G^. 

Prom this theorem we might expect that all asymmetric perfect Delaunay 
polytopes are obtained from symmetric ones by sectioning, and at our current 
state of knowledge this appears to be the case. This would be a fortunate 
turn of events, since the growth of classes of symmeteric perfect Delaunay 
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d 


2 (("..) + (,.:,)) = 2(::,') 


centrally-symmetric 




d-1 


d(d+l) ^ 
2 


asymmetric 



Table 1: Properties of constructed perfect Delaunay polytopes 



polytopes with dimension is slower than for all other geometric phenomena 
associated with point lattices that we know about. 

We summarize properties of the constructed polytopes in the following 
table (in this notation groups of entries separated by semicolumns can only 
be permuted between themselves). 

The following table gives coordinates of the vertices of G-topes, discovered 
in 2001 by Erdahl and Rybnikov. 



[O'^] X 1 




-l,0''-i;-l] X (rf- 1) 




l'^-^-(rf-3)] X (rf- 1) 


[0.r'-2: (,/_,!)] X (d-l) 




[i.(y'-ho\ X {d- 1) 



Polytope Pj]^ is afiinely equivalent to Gosset's 821 = G . Polytope is 
affinely equivalent to Gosset's 221 = G^. 



4 Delaunay Property Part of Main Theorem 
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We will use the following notation: 0^ = {Ylt=i > ^2(^) 
02 (x) is the squared Euchdean distance from x to the hne hn j. The following 
theorem proves the Delaunay property of polytopes P^f, asserted in Theorem 
4. 

Theorem 6 Lets>l,k>2, andd> A;(2s+1)+1. Letf = [{-lf''^,l<^-^/% 
A = < ei, . . . , Cd, A° = {A e A : A • Z° = 1( (mod 2))}. Then there is 

a positive definite quadratic form of the type 

0(x) =a0i(x)+/?02(x), (1) 

were a, /5 > 0, such that the polytope Pf j^ is a Delaunay polytope in the affine 
lattice A" with respect to 4>. 
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Here n — d — k and ei,. . . ,6^ stand for the canonical basis of R*^. 



Lemma 4.1 Suppose 0(x) = q;0i(x) + /302(x) where a,f3 > 0. Then all 
points A e A° which are closest to with respect to (j), ie., 

0(A) = imn{(f){u) : u G A°} 

are, up to permutations of components, of the type [l', O*^"'''] + a^, where 
— I < Z < |; a e Z. Here l' /or a negative I means I times —1. This 
representation is unique. 

Proof. Suppose A e A° is closest to with respect to in A°. Let A = 
aiCi + • • • + a^ed + a-, ai, . . . , aa, a & 'Z, m — ai + • ■ ■ + Ud. We have 




d 
i=l 



2 J 

» 2 

m 



The numbers are of two consecutive integer values. Otherwise there 
would be Oj and aj such that Oj — Oj > 2. Consider the vector A' — aiei + 
• • • + (oj — l)ej + ■ — h (flj + l)ej + ■ ■ ■ + ttdCd + a^. We have 02('^') ~ ^2(^) = 
(oi - 1)2 + (a^- + 1)2 - - a^. = 2{aj - a^) + 2 < -2 and j ■ A' = j ■ A. Since 
A' G A° and 0(A') < 0(A), it follows that the vector A is not closest to 
which is a contradiction. 

Now let b be the smallest of the values of numbers Oj. Subtract b{ei + 
■ ■ • + ed) from the first part and add an equal value of bk^ to the second part 
of the existing representation of A. After a permutation of components, the 
vector A is equal to [I'-C^"'] + (a + bk)j-^ where < I < d. If / > |, again 
subtract [I"'] from the first summand and add j to the second summand to 
get the required representation. 

To prove uniqueness, note that the components Aj of the vector [l', 0''"'''] + 

are of at most two values. In vector [l', C^^'''], either Is fill the positions 
of the bigger Aj, or —Is fill the positions of the smaller A^. Since — | < / < |, 
the choice is unique. □ □ 



Lemma 4.2 All vectors in the affine lattice A° which are closest to with 
respect to a form = a0i(x) + /302(x), a,P > 0, belong to the set {A e 
A° : |A ■ j| < -}. If X is closest to with respect to and |A • j| = -, then 
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Proof. By multiplying both a and f3 by the same positive number, we can 
assume that the minimal value of d) on A'' is 1. Since - G A°, (hi-) = a(-Y > 
1. Let A G M'^ be a point with 0(A) = 1. Represent A = ^ j + m where 7 G M, 
M . j = 0. We have 0(A) = + = 1, therefore < 1. Since 

c^iiY > 1) we have proven that I7I < 1, so |A ■ j| = \ < ^. 

If the latter inequality holds strictly, then I7I = 1 and 0(A) = + 
= 1. Since a(^)^ > 1, we necessarily have /3 = so A = □ 

Proof of Theorem 6: 

We define a collection of points in A'' which contains, up to sign and 
permutation of components, all vectors of the affine lattice A" which have 
the minimal distance to with respect to any quadratic form of the type 
(P). Then we pick the form so that the vertices of P^^ are minimal vectors 
and other points of M are not. Below is an implementation of this program. 

(A) Consider the set 

M = (a = [1', 0"-l'l] + a i G A° : < A ■ j < -, -- < / < -1 U {^}. (3) 

n n 2 2 ) n 

By the two previous lemmas, for every minimal vector A G A° of the quadratic 
form 0, either A or —A, after permutation the components if necessary, be- 
longs to this set. 

For A = [1', 0'^"'] + a^, we have 

A.j = /+ — ,0,(A)= + — ^ ,<P^{X) = \l\-^^,f-X = a+l (mod 2). (4) 
n \ n J a 

From the calculations above, we get 

M = I [1', O'^-I'l] + : / + a = 1(2), < /A; + arf < rf, -- < / < - 1 U {^}. 

(5) 

We define a mapping from M.'^ to by (pi 2(x) = (0]^(x), 02 (x)). We will 
call the image of M the diagram. Lines parallel to 0j^-axis in will be called 
horizontal. 

(B) For Ai, A2 G M, points (pi 2(^1)5 0i 2(^2) belong to the same horizontal 
line if and only if Ai = ±A2. If A and —A both belong to M, then A ■ j = 0. 



15 



To prove that, take A = [1^ 0<^-l'l] + G M. If / 7^ 0, then the condition 

< Ik + ad < d uniquely defines a from a given value of /, and if Z = 0, then 
a = 1. Therefore / uniquely defines a. 

The function Z — >■ |Z| — ^ is even and increasing on [0, |] so two different 
points of the diagram may belong to the same horizontal line if and only 
if their preimages are Ai = [l',0'^~'''] + ai- and A2 = [1^', O*^^''!] + 02- for 
some l,ai,a2. If Z 7^ 0, then < Ik + aid < d, < —Ik + a2d < d. Adding 
these inequalities, we get < (ai + a2)d < 2d, so ai + 02 G {0, 1}. If 

01 + 02 = 0, then Ai = — A2 so (p^ 2(^1) = 0i 2(-^2)- If cti + 02 = 1, then the 
numbers / + ai and — / + 02 have different parity which contradicts conditions 
/ + ai = — / + 02 = 1(2). If / = 0, then Oi = 02 = 1. This proves our claim. 

(C) Our method of constructing perfect Delaunay polytopes can be sum- 
marised as follows. We note that each line axi + [3x2 = 1, where > 0, 
which supports the edge of the convex hull of the diagram, gives rise to a 
quadratic form 0(x) = a0]^(x) + /302(x) such that the ellipsoid 0(x) = 1 
circumscribes the points of the affine lattice hP which fall into the endpoints 
of the edge and has no lattice points inside. 

(D) First consider case ci = 7, /c = 4, s = 1. The diagram for these values 
of parameters is shown in figure El We see that the line |xi + |x2 = 1 passes 

through points ^i^aKi) = 0i,2([l>O6]) and 01,2^,2) = 0l,2([-l^O5] + i), 
and all other points of the diagram are contained in the open half-plane 
+ 1x2 > 1. This means that polytope Pj^ is a Delaunay polytope with 
respect to quadratic form |0i(x) + |02(^) ~ 1^ ' ^ + |(j ' ^)^- 

(E) Now suppose that d>S,A<k< rf/2, k = 0(2) and 1 < s < f - 1. 
Suppose A = [l',0'^-'] + ai e M and < 02(A) < i - (or> equivalently, 
|/| < f). If |/| < f , then / > and 

A= [l',0^-']-(/-l)^. (6) 
If |/| = I (this implies that | is an integer number), then 

±x = [ii,Q^-i]-d-i)l (7) 

k n 

To prove this, first check that these points belong to set M, and then use 
statement in paragraph (B) that there is at most two points A = [l', C^^'''] + 
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G M with a given value of |/| , and there are two if and only if Ik + ad = 0, 
in which case the points are ±A. 

We'll use notation v^^ = 0"^^'] - (/ - 1)^ for < s < f . We have 




All points 01^2 ("^fc.s) therefore belong to the curve 

(9) 

which touches the vertical axis in the point (0, f - ^) when t = f . The curve 
is drawn in dashed line in figure 01 The portion of the curve for < t < | is 
the graph of a convex function. 

Therefore, for each 0<s<-^ — Iwe can find a line with equation 
axi + [3x2 = 1 which passes through points (t>i2iyks)^ '/'i 2(^^5+1)5 supports 
the convex hull of the diagram and does not contain points 0^ 2('^) ^ ^1 
A 7^ i'Wfc itif Quadratic form 0(x) = a0]^(x) + /302(x) defines an empty 
ellipsoid with center which contains vertices of P^^ on its boundary and 
does not contain any other lattice points. 

We have proven that ^ is a Delaunay polytope in affine lattice Aq with 
respect to a unique quadratic form </> = for d > 7. Explicit formula (jH) 
for 0^ g is established by a trivial calculation. 

□ 

Example of the diagram for d = 19, A; = 6 is shown in figure 0] (right-hand 
image). Note that not all points belong to the curve ©. 

5 Perfection Property Part of Main Theorem 

Theorem 7 Let d, k and s be as in Theorem ^ Then, there can he at 
most one positive definite quadratic form whose ellipsoid 0(x — c) = 1 
circumscribes the polytope P^ ^ for some c G M'^. 

Proof. As before, n = d — k. Since for any given positive definite quadratic 
form the center of the its ellipsoid which circumscribes P^^ is defined uniquely. 
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and since the polytope P^^ has as the center of symmetry, is the cen- 
ter of the circumscribing eUipsoid. Therefore, if (f) and ip are two positive 
definite quadratic forms which circumscribe P^^ by elhpsoids of radius 1, 
then the elhpsoids are given by equations 0(x) = 1, ijj{'x) = 1. Hence 
(0 — ip)\vd ^ = 0. We consider an arbitrary quadratic form / such 
that flvd = and prove that f — 0. 

k,s k,s+l 

We will use the following symmetrization technique. Suppose that G is 
a subgroup of the group Sd of permutations on d elements Id — {1, . . . , d}. 
The symmetrization F of the form / by G is defined as 

F{xi, ...,Xd) = ^ f{Xal: ■ ■ ■ , X^d)- (10) 
aeG 

Since polytope P^^ is invariant under transformations of R.^ defined by per- 
mutations of coordinates, all components in the above sum are on the set 
of vertices of Therefore F|vd , =0. 

'^'■5 ' ^k,s'~'^k,s+l 

Firstly we prove that if a form / satisfies the condition flyd yyd = 0, 

k.s k,s-\-l 

then the sum of the diagonal elements of /, and the sum of off-diagonal 
elements are 0. We symmetrize / by the group Sd and get a form F with 
diagonal coefficients proportional to the sum of diagonal coefficients of /, and 
non-diagonal coefficients proportional to the sum of non-diagonal coefficients 
of /. Therefore F can be written as P(x) = q;0i(x) -|- /902(x). Suppose 
that a and (3 are not both equal to 0. Since Flyd ^jyd = 0, the line axi + 

k.s k,s-\-l 

13x2 = passes through the points (pi^iYI^,.) and <Pijyks+i)^ where (j)-^ gl^) = 
(0]^(x), 02 (x)). This is impossible because the line that passes through these 
points is uniquely defined and does not contain 0. This contradiction proves 
our claim. 

Next we prove that all diagonal elements of the form / are equal to 0. 
Suppose that one of them is nonzero. Without a limitation of generality we 
may assume that t = /n 7^ 0. Let F be the symmetrization of / by the 
group G = St{l) permutations which leave 1 fixed. The matrix of F is 



(11) 



t 




/9 




■ P 




a 


S 




. s 


/? 


6 


a 


6 


. 6 




5 




S 


a S 




S 






S a 
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Consider the following vectors: Vi — [1*, O'' — (s — 1)^ e V^^, V2 — 
0'^-*-^] - si e V^^^_^^. The conditions F{vi) = 0, F{v2) = yield 

t + 2{s- 1)P +{s- l)a + {s-l){s- 2)5- 

2(s -l)(t + (s + d-2)P + (s- l)a +(s-l){d- 2)5) 

n 

{s - l)\t + 2{d - 1)P + (rf - l)a +{d- l){d - 2)5) _ ^ 



2s(t + (s + d-l)P + sa + s(d-2)5) 

t + 2s(3 + sa + s{s - 1)5 ^ ^ ^— ^ ^+ 

n 

s\t + 2{d - 1)P +{d- l)a +{d-l){d- 2)5) _ ^ 

n? 

We also have the conditions that the sums of the diagonal and the off- 
diagonal elements are equal to 0: 

t+{d-l)a = 0, {d -l){d- 2)5 + 2{d - l)(3 = 0. (14) 

The determinant of the system of these 4 equations in variables a, (3, 5 and 
t is equal to 

-{d-l){s-d + l)(s-d)(d-2-n) (15) 
n 

and it is not equal to for the specified values of d, n and s. Hence t = 0. 

Next we prove that all off-diagonal elements of / are equal to 0. Without 
a limitation of generality we can assume that /12 7^ 0. Let F be the sym- 
metrization of / by the group of permutations which map the set {1,2} onto 
itself. The matrix of F is 



t 








■ P 


p 


a 


5 




. 5 


p 


5 


a 


5 


. 5 




5 




5 


a 5 




5 






5 a 



(16) 
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where a 7^ 0. We consider the following vectors: Vi = [1*"'"^,0'^ ^ ^] — G 
V;^^^ and V2 = [0'^-'-\ I'+V^i e V^^^^. The equations F{vi) = 0, F{v2) = 
yield 



A{s-l)(3 + 2a + {s-l){s-2)5-2-{2a + 2{s-l+d-2)p + {s-l){d-3)6) + 

n 

{s - If {2a + A{d - 2)13 +{d-2){d- 3)6) ^ 



, , 2{s + l)s{2l3 + {d-3)5) s\2a + A{d-2)p + {d-2){d~3)5) ^ 

(s+l)sd \ = 0. 

n 77,^ 

(18) 

We also know that the sum of off-diagonal elements of the matrix of F is 
equal to 0: 

2a + A{d-2)l3+{d-2){d-3)5 = Q (19) 
so, with some simplifications, the previous two equations can be rewritten as 



4(s-l)/3 + 2a + (s-l)(s-2)(5-2-(2a + 2(s + rf-3)/?+(s-l)(rf-3)(5) = 0, 

n 

^ _ 2(2,9 +(d-3M)^^ (20) 

n 

The systems of equations (fTUI) and ipUl) in variables a, /? and 5 has de- 
terminant %{d — 2 — n){—d -|- s -|- 1) which is not equal to for the specified 
parameters d, n and s which proves that the system has only a solution. 
In particular, a = 0. We have proven that all off-diagonal elements of the 
matrix of form / are equal to 0. □ 



6 Proof of Theorem 5 

Assume that D C M"'+\ with ^ aff Consider the lattice A =< vert D >z. 
Let V G vertD be any vertex. Consider polytope D' = v — D. We have 
vertD' C A. One can see that there is a positive definite quadratic form 
whose ellipsoid circumscribes vert D fl vert D' and does not have any other 
lattice points inside or on the boundary. Note that since D satisfies the 
perfection property, is defined up to a parameter and a scale factor. 
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Wc can now "stretch" the elhpsoid until it touches some other lattice 
point (s) X. If the vertices of D cannot be embedded into the set of vertices 
of a centrally symmetric perfect Delaunay polytope, then the ellipsoid will 
indeed meet points in lattice A. Proving that means proving that the cylinder 
which circumscribes both polytopes D and D' contains other points of A. If 
that is not true, project A onto linD' along t. The image Ai is either a 
lattice, in which case the ellipsoid which circumscribes D' will circumscribe 
a centrally symmetric perfect Delaunay polytope in Ai, or the direct sum of 
a lattice and a linear subspace, in which case the set of vertices of D' has 
dimension less than d-a contradiction. 

The resulting ellipsoid circumscribes the polytope D" — conv(vert D U 
vert D'UX) and has no other lattice points inside or on the boundary. Hence 
D" satisfies the Delaunay property. It is easy to see that it also satisfies the 
perfectness property. 

7 Summary 

This paper has presented two equivalent definitions of perfect Delaunay poly- 
tope, one related to the stTidy of edge forms in L-type domains, the other to 
a generalization of perfect forms to inhomogeneous case. Edge forms classi- 
fication is one of the cornerstone problems in geometry of positive quadratic 
forms; as Voronoi noticed in 1908 , the study of minima of inhomogeneous 
forms is strongly related to the study of perfect homogeneous quadratic 
forms. 

We observed that each centrally perfect Delaunay polytope can be embed- 
ded into a centrally symmetric one. Therefore centrally symmetric perfect 
Delaunay polytopes play a special role. We conjectured that the growth of 
the number of perfect Delaunay polytopes in dimension d up to affine equiv- 
alence is not prohibitive; M. Dutour (2004) discovered that there is only one 
perfect Delaunay polytope in dimension 6. We have given a new construction 
of infinite series of perfect Delaunay polytopes in dimensions d > 6, based 
on the Gosset polytope 821. 

Further study of perfect Delaunay polytopes can have two directions: one 
to estimate the growth of number of types of centrally symmetric Delaunay 
supertopes, the other to find higher-dimensional analogues of the 15, 16 and 
22, 23-dimensional examples. 
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